Spin-transfer pulse switching: Prom the dynamic to the thermally activated regime 
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The effect of thermal fluctuations on spin-transfer switching has been studied for a broad range 
of time scales (sub-ns to seconds) in a model system, a uniaxial thin film nanomagnet. The nano- 
magnet is incorporated into a spin- valve nanopillar, which is subject to spin-polarized current pulses 
of variable amplitude and duration. Two physical regimes are clearly distinguished: a long pulse 
duration regime, in which reversal occurs by spin-transfer assisted thermal activation over an en- 
ergy barrier, and a short time large pulse amplitude regime, in which the switching probability is 
determined by the spin angular momentum in the current pulse. 

PACS numbers: 75.78.Jp, 85.75.-d, 75.75.Jn, 85.75.Bb 



Magnetization dynamics of a nanomagnet in the pres- 
ence of thermal noise is a topic of great fundamental in- 
terest and one of importance to magnetic information 
technologies. Of interest is the probability of a nano- 
magnet at finite temperature to reverse its direction of 
magnetization by thermal activation over an energy bar- 
rier. Spin-transfer torques [H-[l] significantly affect ther- 
mally activated magnetization reversal, either enhancing 
or suppressing the reversal rate depending on the sign of 
the current [a, [f|. Further, thermal fluctuations play a 
central role in spin-transfer device characteristics. A cur- 
rent greater than a critical current is needed to reverse 
the magnetization direction at zero temperature (i.e., in 
the absence of thermal noise) [4|. However, spin-transfer 
driven reversal typically occurs for currents far less than 
this critical current because of thermal fluctuations. 

Analytic models for finite temperature spin-torque dy- 
namic have only been developed for uniaxial single do- 
main nanomagnets, i.e. a nanomagnet with two energy 
minima, magnetization "up" and "down", separated by 
an anisotropy barrier 0-(3- However, experiments have 
mainly explored soft thin film magnetic elements with a 
biaxial anisotropy, a shape anisotropy that leads to in- 
plane magnetization with a preferred axis in the plane 
|a-flOl|. The recent development of advanced nanopil- 
lar devices based on perpendicular magnetic anisotropy 
(PMA) materials [ll|, [l2j now permit a direct compari- 
son to analytic models. Further, most studies, including 
our own initial studies of pulse magnetization reversal 
[l3j . examine the mean switching time, rather than the 
full probability distribution. Here we present results for 
spin-transfer switching of a uniaxial nanomagnet subject 
to current pulses from second to sub-nanosecond dura- 
tion and measure the switching probability versus pulse 
duration. 

Experiments were conducted at room temperature on 
spin- valve nanopillars that consist of two magnetic layers, 
a reference and a free layer, both with a strong uniaxial 
anisotropy perpendicular to the plane of the layers. The 
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FIG. 1: (a) Rac vs. Ida of a 100 nm device at room tempera- 
ture in zero applied field, (b) Schematic of pulse experiments: 
a current pulse is applied to the nanopillar either via a bias 
tee or switch, (c) Pulse switching results for pulses of du- 
ration 0.3 ns to 1 s: pulse amplitude for P — 0.5 vs. pulse 
duration. Inset: Switching rate (1/t) vs. pulse amplitude for 
short pulses showing 1/t cc I. The intercept (1/t = 0) is 
I c o- The dashed line in the main figure shows this behavior 
on a logarithmic scale. The dashed-dotted line is a fit to the 
long-time data using Eq. [T] 



free layer is a Co/Ni multilayer with a room tempera- 
ture coercive field of about 0.1 T and the reference layer 
consists of Co/Pt and Co/Ni multilayers with a much 
larger coercive field (~ 1 T). The layer stack [13( is pat- 
terned into 50 and 100 nm square nanopillars with top 
and bottom electrical contacts such that current flows 
perpendicular to the plane of the layers. 22 junctions 
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were studied and here we focus on the characteristics of 
two representative samples. 

Figure 0Ja) shows a measurement of the differential 
resistance versus dc current of a 100 nm device in zero 
applied field. The high resistance branch corresponds to 
a state in which the free layer and reference layer mag- 
netizations are antiparallel (AP) and the low resistance 
branch corresponds to the two layers magnetized in a par- 
allel (P) configuration (MR=0.3%). A current of 4 mA 
leads to single step switching of the device state from AP 
to P. Conversely, —7 mA changes the device state from 
P to AP. 

To determine the switching probability as a function of 
time we apply current pulses of variable amplitude and 
duration, as illustrated in Fig. QJb). The sample state 
was determined by measuring the device resistance using 
a small (< 300 fj,A) ac current and a lock-in amplifier 
before and after the pulse. The same pulse amplitude and 
duration is repeated 100 to 10, 000 times to determine the 
switching probability. Since the state of the free layer 
is very stable in the absence of an applied field and dc 
current, switching only occurs through the application of 
a current pulse. We focus on AP to P switching events 
and have confirmed that similar behavior is observed for 
switching from the P to AP states. 

Figure [He) shows the switching behavior from the AP 
to the P state over nine order of magnitude in current 
pulse duration. The pulse amplitude I corresponding 
to the observation of switching for half of the events 
(P = 0.5) is plotted versus pulse duration r on a logarith- 
mic scale. For short pulse durations (r < 10 ns) the pulse 
amplitude required to switch the nanomagnet increases 
dramatically. For long pulses the pulse amplitude de- 
pends weakly on pulse duration, varying logarithmically 
(r>10/xs). 

The short and long time switching characteristics are 
thus quite distinct. The inset of Fig. [ljc) shows 1/r 
plotted versus pulse amplitude at short pulse durations. 
This boundary follows the form 1/r — A(I — I c o)i which 
is expected based on the conservation of angular momen- 
tum [J] . The zero temperature critical current I c o in this 
expression reflects the portion of the spin-angular mo- 
mentum that is needed to overcome the magnetization 
damping and the parameter A characterizes the link be- 
tween charge and spin-angular momentum transport. We 
find J c0 = 6.55 mA and A = 2.0x 10 11 C" 1 for this sam- 
ple. The dashed line in the main part of Fig. fjjc) shows 
this same fit to the short time data on a logarithmic pulse 
duration scale, showing that this form describes the data 
up to pulse durations of 5 ns. 

At long times the switching boundary was fit to the 
form ifl: 

I = I c0 (l-^ ln(-r/(r ln(l - P))j , (1) 

where Uo is the energy barrier in the absence of the spin- 
current. This expression fits the data for r > 6 /is well, as 
seen by the dashed-dotted line in Fig. [He). The slope of 
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FIG. 2: Switching probability versus pulse duration in three 
regimes: (a) Short-time limit, / > I c o = 6.55 mA; (b) Inter- 
mediate regime, I — I c q. (c) Long-time limit, I < I c o- The 
solid lines are fits to Eq. [2] and the dashed curves are fits to 

Eq.ru 

the curve gives the ratio £ = Uo/kT = 63 and the inter- 
cept of this curve with I c q gives the Arrhenius pref actor, 
To = 24 ps, which we discuss further below. Figure [ljc) 
also shows that at intermediate time scales, from ~ 10 ns 
to 6 /Lts, the switching boundary follows neither the 
short or long time forms, i.e. this is a crossover regime. 

In order to understand the magnetization switching 
processes in more detail we studied the full switching 
probability distributions in the three regimes, i.e. at 
short, intermediate and long time scale (Fig. It is 
clear that the switching probability depends differently 
on pulse duration in these three cases. For short times 
the curves are sigmoidal and for long times the switch- 
ing probability depends exponentially on pulse duration, 
P = 1 - exp(-r/rA), with t a = t exp(£(l - I /Ida)). A 
fit of this data gives To = 24 ps and £ = 63, consistent 
with the fit to the long time data in Fig. Hie). 

In the short pulse time limit the reversal process can no 
longer be considered thermal activation over an energy 
barrier. In this case, the current amplitude is larger than 
the critical current I c o and spin-torque quickly drives the 
magnetization reversal, starting from a thermally dis- 
tributed initial state. The reversal time depends sensi- 
tively on the initial magnetization state and the distribu- 
tion of initial magnetization states leads to a distribution 
of reversal times. This limit was considered theoretically 
in Ref. Q for a uniaxial single domain nanomagnet. The 
following expressing for the switching probability was de- 
rived: 

P = exp(-4«expplW^))), (2) 

where td = l/( a 7Mo-f^fe) is the characteristic time as- 
sociated with the spin-transfer driven magnetization dy- 
namics, a is the Gilbert damping, 7 is the gyromagnetic 
ratio and Hk is the anisotropy field. This expression fits 
the data well, as seen in Fig. Ufa), with £ = 63 and 
td = 260 ps. 

At intermediate times (Fig.^b)), the switching proba- 
bility has a distinct form, it is neither the sigmoidal form 
measured at short times (the solid curve in Fig. Hfb)) 
nor the exponential form (the dashed curve in Fig.[2][b)) 
observed at longer time scales. On such time scales ther- 
mal fluctuations during magnetization reversal and spin- 
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FIG. 3: Switching probability as a function of the scaled an- 
gular momentum ,At(I — Ico), m the short time regime for 
both 50 and 100 nm devices. Both data taken with variable 
pulse duration and variable pulse amplitude are shown. For 
t < 5 ns P only depends on the net spin angular momentum 
in the pulse. 



torque driven processes are both important-and cannot 
be separated. 

An interesting prediction of the model for short time 
switching (Eq. [2]) is that the switching probability only 
depends on the total spin angular momentum in the 
pulse, phlr/e, where p is the spin-polarization of the 
current Q. We test this model and more generally 
whether the spin angular momentum is the relevant ex- 
perimental variable by plotting the switching probabil- 
ity in Fig. [UJa) versus the scaled angular momentum, 
At (I — I c q) (Fig. [3]). Included in this plot are additional 
measurements of the switching probability in which we 
varied the pulse amplitude for a constant pulse dura- 
tion. Figure [3] shows that there is a good scaling of the 
data, particularly for short pulse durations. Deviations 
from the scaling form (Eq. [5]) occur for pulse durations 
of 5 ns and greater. This indicates the time scales on 
which thermal fluctuations influence the switching pro- 
cess during the current pulse. Pulse studies conducted on 
a 50 nm device are also shown in Fig. [3] For this sample 
J c0 = 1.04 mA and A = 1.5 x 10 12 C _1 . Data in this 
case, with pulse amplitudes / > 2I c q, closely follow the 



scaling curve. 

The single domain model assumes an energy den- 
sity: E = —iioM s Hkm 2 z /2 1 where m z is the normalized 
magnetization in the z-direction. The energy barrier is 
Uq = fioM s HkV/2. fioHk of our nanomagnet is 0.25 T 
giving an energy barrier of £ = 360 at room temperature, 
six times larger than the energy barrier found from the 
data in Fig. [TJc) and Fig. Ufa), (c). This implies that the 
nanomagnet is not reversing as a single magnetic domain. 
Further, the attempt time in a single domain model is 
given by r = 0F(1 + a 2 )/{2 1 ^H k a^l) 0,111. With 
£ = 63 and a = 0.01 we find tq = 1.6 ns and with 
a = 0.1, To = 162 ps. In either case, this time is sig- 
nificantly larger than the 24 ps determined in our exper- 
iment. While the To we find is small compared to the 
macrospin model, it is within the range of values deter- 
mined in other experimental studies. Krause et al. (l6j 
report even shorter tq, from 10 -16 to 10 -13 s, for nanos- 
tructures studied by scanning tunneling microscopy. In 
the short time limit the single domain model predicts 
that td = {ajpoHk)^ 1 ■ Taking a — 0.01 to 0.1 gives 
Td = 1.4 to 14 ns, far larger than td deduced from the 
data in Fig. [D(a), td = 0.26 ns. These comparisons show 
that the microscopic time scales are shorter than those 
predicted in the single domain model. 

Our results demonstrate that there are two distinct 
spin torque switching processes, spin-transfer assisted 
thermal activation at long times and one dominated by 
angular momentum conservation at short times. The 
measured switching probability distributions at both 
short and long times are in good agreement with a sin- 
gle domain model that includes spin-transfer torques and 
thermal fluctuations 043 • However, parameters such 
as the energy barrier deduced from fits to these models 
are far less than that expected. This suggests that the 
switching occurs by rotation of small part of the nano- 
magnet, whose dynamics is nonetheless captured by a 
macrospin model, as was found in field driven reversal of 
nanostructures with PMA [l7l ]. 
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